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Abstract 

The energy spectrum of primary cosmic rays is explained by 
particles emitted during a thermal expansion of explosive objects Inside 
and near the galaxy, remnants of which may be supernova and/or active 
galaxies, or even stars or galaxies that disappeared from our sight 
after the explosion. A power law energy spectrum for cosmic rays, 

E" 01-1 , is obtained from an expansion rate t « R a . Using the solution of 
the Einstein equation, we obtain a spectrum which agrees very well with 
experimental data. The implication of an Inflationary early universe on 
the cosmic ray spectrum is also discussed. It is also suggested that 
the conflict between this model and the singularity theorem in classical 
general relativity may be eliminated by quantum effects. 
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Since the discovery of cosmic rays early this century, an impressive 
amount of experimental data has been accumulated. Yet, the origin of cosmi< 
rays defies the understanding of physicists. Several important questions 
are: What is the fraction of galactic and extragalactic components of 

primary cosmic rays? How can one understand the power law energy spectrum 
(~E-2*5 - E“ 3 for the energy range 10 9 eV < E < 10 20 eV) ? How do they 
attain such high energies? Despite various attempts in the past to answer 
these questions, we are still left in the dark. 1 In this article, we 

present a model which explains some of the features of cosmic rays describe* 
above. 

The following are the assumptions which we propose to make: 

(a) The sources of primary cosmic rays are supernova, galaxies and/or 
clusters of galaxies which have exploded in the past, after the 
gravitational contraction. Active galaxies such as radio- or X-ray 
galaxies, Seifert galaxies, or QSO's may be the remnants of such an 
explosion. Or some of them might have disappeared completely from our 
sight after the explosion. 

(b) The explosion of such an object is a replica of the expansion of the 
universe at a smaller scale. It goes from expansion in a radiation 
dominated era at extremely high temperature to that in a matter 
dominated era at lower temperature. It may even undergo a phase like 
that of an inflationary early universe. 

(c) In the course of expansion, high energy particles of various kinds are 
emitted because of nonequilibrium processes. The energy distribution 
of the emitted particles reflects the black body radiation law at the 
temperature T at which the particles are emitted. The number of such 
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particles is not large and the system is approximately in equilibrium 
at each instant since the collision time is much smaller than the 
expansion time. 


Under the assumptions, the number of particles of type x emitted with 
energy E is given by 


f ( E) -ils+il / n (E/kT) 


2u 


E 2 v s dt 

E/kT - p/kT 


± 1 


( 1 ) 


where V$ is the effective volume around the surface of the system which 
emits the particles, n x (E/kT)dt is the fraction of particles x emitted in 
time interval dt and p and s are the chemical potential (zero) and the spin 
for fermions (bosons) of type x. The +(-) sign in the denominator is for 
fermions (bosons). The volume V$ is taken to be 


V$ - 4irR^d = 4 tt 


a? 

(kT) 3 


( 2 ) 


since the surface of the particle distribution is defined with uncertainty 
ax = d = 1/kT, where a is a constant defined by R = a/kT. For an expanding 
system 

t « bR“ , (3) 


where b is constant and 2 


<x=2 for the radiation dominated era 


(4a) 


and 


o=3/2 for the matter dominated era (4b) 

(the Ei nstei n-deSitter flat space) 

The function n x ( E/kT) is unknown, but is assumed to scale as a function of 
E/kT. The result is essentially the same if the E (or T) dependence of 
n x (E/kT, E (or T) ) is slowly varying except in E/kT. The chemical potential 


for fermions is obtained by the condition 
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or equivalently 


N = (2s+l) 


2 it 


(kT)* 


J 

o 


x 2 dx 


g x - u/kT + j 


u 0 = u/kT = 


9 (~ 


V(kT) 


2*2 

3 (2s+l) 


)» 


(5) 


( 6 ) 


where N is the total number of particles in the system and V is the volume 
given by 


V = i* r3 = il a3 
3 3 (kT ) 3 


(7) 


(At high temperature, g(x) = £,n(x/ 2 ).) Obviously, y 0 is independent of 
temperature since VT 3 is constant during the course of the expansion. (u o =0 
for bosons, as was stated earlier.) Using eqs.(l)-(4), we obtain 


fx(E> 


X, a 

-o+l 


( 8 ) 


where 


. _ 2 ( 2 s+l)gb(a) 2+ a 7 n x (s)s dS 

* 6 e S_u ° ± 1 


x,a 


(9) 


is a constant. Then, the total energy spectrum is given by 


f (E) * 


frg+T 


where 


k - x 

X 


( 10 ) 


(ID 


The power law energy spectrum ( 8 ) and (10) with the values of a in eq.(4). 


is 



